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that the classical nature of the flow equation allows to eliminate all cutoff effects at 0{a?) 
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modified initial condition for the flow equation. A consistency test in perturbation theory 
is passed and allows to determine all counterterm coefficients to lowest non-trivial order in 
the coupling. 
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1 Introduction 


In recent years the Yang-Mills gradient flow has been established as a very promising 
new tool to study non-perturbative aspects of strongly coupled gauge theories [1-4]. The 
gradient flow defines a deterministic mapping from the original gauge field A^{x) to a 
smoothed gauge field configuration, at flow time t, which is obtained as the 

solution of the gradient flow equation (see appendix A for a summary of our conventions), 

x) = ^ DyGy^{t, x), S^(0, x) = A^(x), (1.1) 

V 

where + [Bf^, •] denotes the gauge covariant derivative and 

Gfj^u = d^Bi^ — dyB^ + [B^i By\^ ( 1 - 2 ) 

is the associated field strength tensor. The name relates to the fact that the right hand 
side of (1.1) is equal to minus the gradient of the Yang-Mills gauge action. Hence, with 
increasing flow time t, the solution, B^{t,x), is driven towards a minimum of the action 
and thus approaches a smooth classical field configuration. 

There is quite some freedom when translating the gradient flow equation to a Euclidean 
space-time lattice. A simple possibility is to choose Wilson’s plaquette action, 5w, and to 
define the lattice gauge held at hnite how time, x), as the solution of the Wilson how 
equation, 

[dtVf,{t,x)] V)j(t,x)^ = -gQ^x,^S^^l[V], (1.3) 

where denotes the Lie-algebra valued derivative with respect to Vf^{t,x). It should be 
noted that similar smoothing operations have long been successfully applied in lattice QCD. 
For example, the stout link smearing technique of ref. [5] can be understood as discretized 
how time version of Eq. (1.3), The essential new element is a theoretical understanding of 
the renormalization properties of the Yang-Mills gradient how. In particular, in [3, 4] it was 
proved to all orders of perturbation theory that QCD at hnite how time t is renormalized 
once it is renormalized at how time t = 0 through the usual renormalizations of the gauge 
coupling and the quark mass parameters. Furthermore, gauge invariant helds at positive 
how time are automatically renormalized and do not mix with other helds of the same 
or lower dimensions. These properties allow to dehne a new class of renormalized gauge 
invariant observables which can be used to probe the theory in various ways. It also opens 
new ways to dehne renormalized composite operators at zero how time; the study of Ward 
identities at positive how times [4, 6, 7] and the applications of the so called “small how 
time expansion” have received much attention recently in this context [8-10]. 

Many current lattice QCD applications of the gradient how only involve the simplest 
possible gauge invariant held, the action density, 

E{t, x) = tT{G^y{t, x)G^y{t, x)}. (1.4) 

fj,,u 

As initially proposed in [2], the expectation value {E(t, x)) can be used for a non-perturbative 
dehnition of either a reference scale or a coupling constant. This has proven very attractive: 
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in large volume simulations it leads to the most precise determination of a reference scale 
(for a recent review cf. [11])- On the other hand, when considered in a finite space-time 
volume the scale evolution of the corresponding coupling [12-16] can be traced with high 
statistical precision (see [17] for a recent review). 

Notwithstanding these nice properties a major practical problem are the relatively 
large cutoff effects which have been observed in several applications (cf. [17] and references 
therein). On general grounds, the leading effects are expected to be of order a?. Their size 
depends on the detailed choices made when translating the flow equation (1.1) to the lattice, 
but also on the discretization of the observable and on the lattice action. Alternative flow 
equations have been tried, e.g. in ref. [18] where the Wilson action was replaced by the 
tree-level improved Liischer-Weisz action, 5 lw [19) 20]. For some attempts to reduce cutoff 
effects in the particular observable {E{t, x)) cf. refs. [21, 22]. Here we would like to proceed 
more systematically by applying the Symanzik procedure [19, 23] to the 4-|-1-dimensional 
local formulation of the theory [3, 24]. This will lead us to a particular choice for the lattice 
flow equation, referred to as the “Zeuthen flow” and defined by 

{dtV^{t,x)) F^(t,x)t = -gl (^1 + , V^{Q,x) = U^,{x ). (1.5) 

Here and V* are the lattice forward and backward covariant derivatives, respectively. 
We will show that the integration of the Zeuthen flow equation does not generate any cutoff 
effects at O(a^). If combined with classical O(a^) improvement of the observable all O(a^) 
effects are eliminated apart from those corresponding to local counterterms in the action 
at zero flow time. We will give a complete list of such counterterms and test our framework 
to lowest non-trivial order in perturbation theory. 

The paper is organized as follows: In Section 2 we recall the definition of the 4-|-l- 
dimensional local theory, with flow time as the added dimension. In Section 3 we discuss 
the general Symanzik procedure and the simplifications due to the special properties of 
this theory. We present the classical a-expansion of both the flow action and the gradient 
flow observable E(t,x), as part of the simplified Symanzik procedure, and carry out the 
standard Symanzik analysis for the O(a^) counterterms at the t = 0 boundary. Section 4 
presents a number of perturbative tests of the O(a^) improved theory, and Section 5 our 
conclusions. We have included three appendices regarding our notations and conventions 
(Appendix A), some technical details pertaining to the classical a-expansion (Appendix B), 
and some explicit expressions used in Section 4 (Appendix C), respectively. 

2 Lattice gauge theory in 4+1 dimensions 

The gradient flow equation can be viewed as a way to define a particular class of observables, 
i.e. fields which are functionals of the fundamental gauge field U^{x). The flow time 
thus appears as an additional parameter which measures the range in space-time over 
which the fundamental gauge field enters into an observable defined in terms of the flowed 
gauge field V^{t,x). The flow time t has dimension length squared and the “smearing 
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(a) So 



Figure 1: Wilson loops of type Sq,Si,S 2 and 53 . 


radius” rt = is usually taken as the corresponding length scale^. Thus, gradient 
flow observables are non-local objects from the perspective of the 4-dimensional gauge 
theory and their properties under renormalization are difficult to assess. Moreover, the 
non-locality prevents a straightforward application of the Symanzik expansion, which is 
our main theoretical tool for understanding the cutoff dependence of the theory. For this 
purpose, it is therefore highly beneficial to follow [4] and view the theory from a 4-|-l- 
dimensional perspective, with flow time as the added dimension. In this re-formulation 
locality is restored in the 4-|-1-dimensional sense, and dimensional counting can be applied 
to classify counterterms to the action and observables. 

We start with the formulation of the lattice set-up, including the introduction of a 
flow-time lattice. The latter should be regarded as an intermediate regularization which 
helps to resolve certain technical issues [4] . While none of this is original it serves for later 
reference and to hx our notation. 

2.1 The 4-dimensional lattice action 

On-shell O(a^) improvement of the 4-dimensional gauge theory can be achieved by intro¬ 
ducing, besides the 4-link plaquette action, further 6 -hnk Wilson loops with appropriately 
chosen coefficients [19]. We will consider a general class of lattice gauge actions parame¬ 
terized by the coefficients Ci{i = 0,1, 2, 3), defined by, 

1 ^ 

= E (2.1) 

^0 j=o we5i 

where the second sum extends over all oriented Wilson loops of type Si. As illustrated in 
Fig. 1, these Wilson loops are the usual plaquettes, Sq, the 2x1 planar loops or “rectangles”, 
5i, the bent rectangles or “chairs”, S 2 , and finally the “parallelograms”, ^ 3 . 

It is conventional to normalize the action by requiring 

Co + 8 ci -|- 16 c2 -|- 8 c 3 = 1 , ( 2 . 2 ) 

such that the standard continuum Yang-Mills action is obtained in the classical continuum 
limit, with any choice of the 3 free parameters^. Popular choices are the Wilson plaquette 
(W) action (cq = 1, 01 , 2,3 = 0) and the tree-level improved Liischer-Weisz (LW) action 
(co = 5/3. Cl = -1/12, C2,3 = 0). 

^The radius rt = amounts to 2 standard deviations in the Gaussian smearing function which appears 
in the relation between B^(t,x) and Afj_{x) to leading order in the coupling. 

^Note, however, that the choice of the coefficients is not completely free, with some constraints arising 
from positivity [19]. Our conventions differ from this reference by the exchange C2 C3. 
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2.2 4+1-dimensional set-up 


Given the 4-dimensional action, the flow equation is now incorporated in the action as a 
constraint, by introducing the Lagrange multiplier field L^(t,x), which is hermitian and 
such that is Lie-algebra valued. The 4+1-dimensional action of this theory then 

takes the form 

roo 

5[y, L] = S^[U, {q}] - 2a^ / dt V tr{L^(t, x)F^{t, x)}, (2.3) 


where the boundary condition. 


1/^(0, x) = Uf,{x), 


(2.4) 


is assumed and 


Ff,{t,x) = a ^ {dtVfj,{t,x))Vf,{t,x)^ - a ^g^dr^^f^SglV], 


(2.5) 


is a shorthand notation which allows to write the lattice gradient flow equation in the form 
Ff^(t,x) = 0. The action Sg[V] is some 4-dimensional lattice gauge action for the flowed 
field Vfj,{t,x), the Wilson action being the simplest choice [cf. eq. (1.3)]. In any case it is 
unrelated to the gauge action Sg[U] in (2.3). How to best define Ff^{t,x) is at the core of 
this work and will be discussed in the next section. 

Given the action for the 4+1-dimensional half space t > 0, expectation values of 
composite fields 0[V, L] are defined as usual, 

(0)=Z-iy D[V]D[L]0[V,L]e^pi-S[V,L]), (1) = 1. (2.6) 


A few remarks are in order: first, the integration over the gauge field V^{t,x) includes the 
integration over its boundary values at t = 0 i.e. the standard 4-dimensional gauge field 
Ufj_{x). Hence, for observables which only depend on [7^, the functional integrals over V^\t>o 
and Lf^ cancel between numerator and denominator, reproducing the standard expectation 
value of the 4-dimensional theory. To see this more explicitly it is convenient to pass to a 
flow time lattice with spacing e and lattice points t = ne [4], 



dt ^ tr{L^(t, x)F^{t, x)} 

X,(l 


37)+/, (t, x; e)}, 

n>0 


(2.7) 


where we have assumed the discretization, 

asF^{t,x;e) = + e,x)Vi,{t,x)^ - eyip (^-gQ-^dx,^,Sg[V]j , (2.8) 

with the correct e —)• 0 limit. Inserting this representation of the action into the functional 
integral, the integration over the fields x) produces a string of +functions^ 

OO 

n n . (2.9) 

x,fi n=0 

These can be eliminated one by one, by integrating over V^{n£,x) for strictly positive n, 
leaving the unconstrained n = 0 integration over the fundamental gauge field intact, as 
expected. 

®For a more careful discussion of the limits involved cf. [4] 
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2.3 Gauge symmetry 

By construction, the 4+1-dimensional action is gauge invariant under t-independent gauge 
transformations, 

x) K{x)V^{t, x)A{x + a/i)'*', (2-10) 

where A(x) is an SU(A^)-valued gauge function. This leads to the transformation, 

x) A{x)F^{t, x)A(x)^ (2.11) 

so that gauge invariance of the action is guaranteed provided that 

Lfj,{t,x) ^ A{x)Lf_i{t,x)A{x)\ (2.12) 

i.e. Lfj_{t,x) must be in the adjoint representation of the gauge group. The integration 
measure is invariant under such a change of variables, so that the gauge symmetry of the 
4-dimensional boundary theory is inherited by the bulk theory. 

It is occasionally useful to generalize the gauge symmetry to the flow-time coordinate 
t, i.e. admit flow-time dependent gauge functions A{t,x). In the continuum theory this 
amounts to replacing t = X 4 , dt —)• D 4 = 84 , + ■] and x) —)• G 4 ^(x 4 , x) [1]. In the 

presence of the lattice cutoff (but continuous t = X 4 ) we define the covariant X 4 -derivative 
by 

V 4 V'^(x) = d 4 Vfj,{x) + B 4 {x)V^{x) - Vfj,{x)B 4 {x + a/i), (2.13) 

where x = (x 4 ,x). This, together with the transformation under an X 4 -dependent gauge 
transformation 

B 4 {x) —> A{x)B 4 {x)A{x)^ + A(x)94A(x)^, (2.14) 

leads to the left hand side of the covariant flow equation transforming as 

[V 4 V;(x)] F^(x)t ^ A{x) [ViV^ix)] y^(x)tA(x)t. (2.15) 

Rendering the t-derivative covariant in the definition of T)j (2.5) one then obtains, 

F^{x) A{x)Ff,{x)A{x)\ (2.16) 

so that 4+1-dimensional gauge invariance is established, provided that L^{x) transforms 
just like F^{x). Discretizing the flow-time coordinate is also straightforward, one just needs 
to elevate the fourth component of the gauge field B 4 to a link field V 4 , with corresponding 
changes in the covariant derivative and gauge tranformation behaviour. 

Finally we note that the Yang-Mills flow equation in the continuum can be written as 

3 

= (2.17) 

u=0 

which shows that the 4+1-dimensional theory, while exactly gauge invariant, does not enjoy 
any generalized Lorentz-symmetry. This is of course already clear from the dimensions, in 
particular, dt and thus B 4 must have mass dimension 2 , in contrast to the usual derivatives 
and gauge helds in 4 dimensions. 
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3 Symanzik improvement to 0{a‘^) 


3.1 Generalities 

The re-formulation of gradient flow observables in terms of a local 4+1-dimensional lattice 
gauge theory creates the standard situation to which Symanzik’s effective theory [23] can 
be applied in the usual way. We start with Symanzik’s effective action which is given as 
an expansion in powers of a^, 

Ses[B,L] = S^°'^^[B,L] + a^S 2 fi[B,L] + a^S 2 AB,L] + 0{a^). (3.1) 

One might worry about odd powers of a arising in a 4+1-dimensional theory. However, as 
we will show in detail in Appendix B, gauge invariance, reflection symmetries and the fact 
that flow time parameter t has mass dimension —2 imply that non-trivial counterterms to 
the action must be even-dimensional. In Eq.(3.1) we have separated the effective action of 
the flow in the 4+1-dimensional volume, from the action 82 ^ with support restricted 
to the 4-dimensional boundary at t = 0. Both parts will be discussed in turn below. Besides 
the effective action, also local observables are described by an effective continuum field, 
again expanded in powers of o^. For a generic local observable O we write 

0es = 0o + a^02 + 0{a^). (3.2) 

To O(a^) the Symanzik expansion of lattice expectation values then takes the form, 

= (Oo) + a^( 02 } - a^(OoS 2 ,{i}c - a^{OoS 2 ,b)c + ^(o^). (3.3) 

Here, the expectation values on the RHS are defined in the continuum theory with respect 
to the continuum action and the notation {■)c serves as a reminder that only the 

connected part contributes to the correlation functions with counterterm insertions, for 
instance 

{OoS 2 ,fi)c = {OoS 2 fi) - (Oo) (52,fl). (3.4) 

As the next step in the Symanzik procedure one determines a basis of counterterms 
both for the action and the observables of interest. In the case of the action these take the 
form 

/'OO r 

S 2 fi[B,B\= / (3.5) 

Jo J 

/ rib 

d^xY,Oi{x), (3.6) 

i=i 

where the fields Qi(t, x) are gauge invariant polynomials in the fundamental fields B^{t, x), 
L^{t,x) and their (space-time and/or flow time) derivatives, and the Oi{x) are similarly 
constructed, but evaluated at t = 0. Since 0 ^ 82 ,^ must be dimensionless the fields Qi must 
have mass dimension 8 and otherwise share all the symmetries with the lattice theory. The 
fields Oi are dimension 6 fields, localized at the t = 0 boundary. One of the important 
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outcomes of the Symanzik analysis are the numbers nfl and Ub of basis elements, where fields 
differing by total (space-time) derivative terms are considered equivalent. Furthermore, 
restricting to on-shell improvement the field equations for and can be used to 

simplify the basis. Given a basis of counterterms the final step of Symanzik’s procedure 
consists in adding lattice representatives of these operators to the lattice action, such that, 
with appropriately chosen coefficients, the terms 52 ,a and S 2 ^b eliminated in Symanzik’s 
effective action for the improved lattice action. 

A similar analysis then needs to be carried out for each observable O of interest, 
i.e. O 2 in Eq. (3.2) is given as a linear combination of local fields of mass dimension 
dim(C)o) + 2 which share all the lattice symmetries with O. While this procedure applies 
to any observables, we will here focus on gradient flow observables, i.e. gauge invariant 
composite fields with support at strictly positive flow times. 

If the full Symanzik procedure as outlined above were really necessary, 0(a‘^) improve¬ 
ment would probably remain an academic curiosity. In particular, a rather long list of 
dimension 8 counterterms for ^ 2,3 could be written down, with little hope for practical 
relevance, so that one might be tempted to give up on systematic 0{a?‘) improvement. 

Before proceeding along these lines, however, it is advisable to have a closer look 
at this particular theory. As shown by Liischer and Weisz, the theory is perturbatively 
renormalizable to all orders in the 4-dimensional gauge coupling g [3]. More precisely, 
if one restricts attention to gauge invariant observables, one just needs to renormalize 
the gauge coupling in the usual way, and also the quark masses if the boundary theory 
is generalized to QCD'^ . Moreover, any composite fields defined at finite flow-time are 
automatically renormalized and do not mix with any other fields of the same or lower 
canonical dimension. The action density (1.4) is a typical example: its renormalization at 
flow time t = 0 requires the subtraction of both a quartic and a logarithmic divergence. 
None of this is required at finite t. It is instructive to consider leading order perturbation 
theory to get a basic understanding of the mechanism at work. Effectively, at finite flow 
time t, integrals over the loop-momentum p are cut off by an exponential suppression factor 
oc exp(—2tp^) in the integrand. This renders most momentum integrals finite, so that one 
is only left with those divergences which are cancelled by the standard counterterms in the 
boundary theory. 

Hence the 4-|-1-dimensional theory enjoys rather special properties. In particular, the 
field plays the role of a Lagrange multiplier field which enforces the gradient flow 
equation as a constraint. The smoothening properties of this equation are related to the 
fact that perturbation theory only generates tree diagrams for the correlation functions of 
gradient flow observables [3]. The Symanzik expansion is then very much simplified as we 
expect the following to hold: 

• The absence of bulk loop diagrams in the perturbative expansion of gradient flow 
observables implies that classical improvement of the flow action yields the O(a^) 
effects exactly, i.e. without any corrections. 

^We assume here that the quark fields only live at t = 0, i.e they are not propagated into the 4-1-1- 
dimensional bulk. For generalizations cf. [4] 
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• By the same argument, non-perturbative 0{a^) improvement of cmoposite operators 
at positive flow time can be achieved by choosing discretizations that do not generate 
0 {a?) effects when expanded classically. 

• The only 0{a?‘) counterterms which receive genuine quantum corrections are the ones 
living in the 4-dimensional boundary at t = 0. The full Symanzik procedure outlined 
above thus needs to be applied only to the t = 0 boundary part, 82 ^ of the Symanzik 
action, and of course to any observable which is at least in part localised at the t = 0 
boundary. 

In the following we hrst remind the reader of the classical a-expansion and then address 
these points in the subsequent subsections one at a time. 


3.2 The classical a-expansion 

According to the preceding discussion the counterterms appearing in S 2 fi and in O 2 for 
gradient flow observables are completely determined by classically expanding the lattice 
action in the 4+1-dimensional volume and the observables under consideration to order a^. 
The classical expansion assumes that the lattice approximates an underlying continuum 
space-time manifold on which a smooth continuum gauge field, x), is defined. The 

lattice gauge field, x), is then related to the continuum gauge field by parallel transport 
along the lattice links. Parameterizing the path along the lattice link from x + aft to x 
by z{u) = X + (1 — u)afL (with parameter u G [0,1]), the precise relation is obtained by 
iteratively solving the differential equation, 

^^ + Bf,{t,z{u))'^v{u) =0, ^(O) = 1. (3.7) 


The solution, v{u = 1) = V)j(t,x), can be concisely written in terms of a path-ordered 
exponential. 


M' 


V^{t,x) = V exp < a / du B^{t, z{u)) 


= l + a du B^{t, z{u)jl) 

Jo 

j-l j-Ul 

+ a^ dui du2 Bfj_{t,z{ui)) Bfj_{t,z{u2)) + 0{a^) 

Jo Jo 

= l + aB^{t,x) + ^a^ {d^Bf,{t,x) + B^(t,x)) +0{a^). 


(3.8) 


(3.9) 

(3.10) 


While it is straightforward to carry out the expansion around a = 0, in practice, even a 
simple gauge invariant quantity like the trace of the plaquette contains 4 link variables 
which need to be expanded and combined to fourth order in a to obtain the leading non¬ 
trivial term. It is therefore highly advisable to perform the expansion efficiently (cf. e.g. [19, 
25]). We here follow Liischer and Weisz [19], who, for fixed indices /i and v, proposed to 
work in the following gauge: 


B^{t,x) = 0 for all x; B^ix) = t) if x^ = 0. (3T1) 
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As a result, the expansion around x = 0 is very much simplified. For example, the plaquette 
field, 

Pfiuit, x) = V^{t, x)Vy{t, X + X + aD) x)^ (3-12) 

is reduced to a single link, 

Piiu{t,0) = Vi^{t,afi) = V exp j duBy + {1 — u)ai'')^ . (3.13) 

Recalling the definition of the path ordered exponential (3.9) one needs the expansion of 
the R-held around a = 0, 

aBy {t, afi + Kav) = a^dnBy{t, 0) + {5^ + iKdfidy] By{t, 0) 

+ + 2K^d^dl] By{t, 0) (3.14) 

+ {<9^ + R^(t, 0) + ..., 

where k is a constant and neglected terms are of order o®. Following [19] the gauge covariant 
expressions can be unambiguously restored, with the result, 

aBy {t, apL + Kav) = a‘^G^y{t, 0) + + InDy} G^y{t, 0) 

+ {Pfi + ^f^DyDfj^ + 3k^Z1^| Gf^y{t, 0) (3.15) 

+—a® {D^ + AnDyD^ + G^y{t^ 0) + ... 

Inserting into the path ordered exponential with appropriate replacements for k, we thus 
obtain the gauge covariant expansion for the plaquette held, 

P^lV = 1 + 0 ?G^y + -a^{D^ + Dy)G^y 

+ + 2F>3) Gij,y + 6 G^yG^y^ 

+ W + 299.1?^ + 2DId^ + Dl} G,y 

+ {(3A>At + ‘2‘Dy){G^y)‘^ + 2G^yDyG^y^ + 0{aP), (3.16) 

which holds for any argument (t,x). Similar expressions can be derived for the other 3 
plaquettes in the fi — u plane: 

Q^y{t, x) = Vy{t, X - oz>) V^(t, X - av)Vy{t, x + ap- aC>)V^{t, x)^ (3-17) 

R^y{t, x) = V^{t, x — ap)^Vy{t, X — ap — ai')^V^{t, x — ap — v)Vy{t, x — az>), (3.18) 

S^y{t, x) = Vy{t, x)V^{t, x - ap + ai')^Vy{t, X - ap)^V^{t, x - ap), (3.19) 

and the next few orders can be obtained with moderate additional effort. 
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3.3 Determination of 52 ,a 

To find the bulk counterterm action 52,a we simply need to apply the classical expansion 
to the bulk action in Eq. (2.3). This essentially amounts to the a-expansion of the gradient 
flow equation, i.e. F^{t,x) in Eq. (2.5). For the first term we find, in the Liischer-Weisz 
gauge (3.11), ^ 

= [ dudtB^{t,{l - u)afl), (3.20) 

Jo 

as all other terms are proportional to = 0. The Taylor expansion can be easily 

performed to all orders in a with the result 


du dtBfj, {t, (1 - u)afi) = ^ 


n=0 


(n + 1)! ^ 


dj:dtB^{t,0). 


(3.21) 


We therefore expect that the correct gauge covariant expression at any lattice point x must 
read 

OO ^ 

a~^ [dtVf,{t, x)] V^{t, x)^ = dtBf,{t, x) + ^ x). (3.22) 

At this point one may wonder whether the gauge covariant expression really follows unam¬ 
biguously from the gauge fixed expansion, in particular, whether the t-derivative always 
has to be to the right of the covariant //-derivatives. That this is indeed correct can be 
established by using the 4-|-1-dimensional gauge symmetry (cf. Subsect 2), which implies 
that the a-expansion of this term must be given as covariant derivatives acting on 

Turning to the second term of (2.5), i.e. the gradient force term, we choose a quite 
general lattice gauge action parameterized by co,i ,2 which includes all 4- and 6-link Wil¬ 
son loops (plaquettes, rectangles, chairs) except the twisted chair s/parallelograms. We 
decompose the action as follows: 


)C0;Cl,C2] — Co5g,pl[E] -|- Cl5g,re[E] -|- C25g,ch[f^]- 


(3.23) 


We first express the gradient force in terms of plaquettes and their covariant derivatives. 
For the plaquette action we then find 


5o^x,/.5’g,pi[E] = ^ (^Pi,„{t,x) -h Q^u{t,xY^ 


AH 


(3.24) 


where we have introduced the projection on the trace-less anti-hermitian part, i.e. for an 
N X N matrix M in colour space we define 


(M)ah = -2tr (r“M) T“ 

For the rectangle action we find, 
go9x,^^Sg^y:e[y] = 

1 / 

+ P^i/(t^ x')S^i/{t^ x) x) 

+ (ttVx)) P^y{t^ x) Qfivit") x)(xV ^u(tf x) 

x)) Q^uit, x) + x)aSIuP^iv{t, x)^ 


(3.25) 


AH 


(3.26) 
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and a similar but slightly more complicated expression is obtained for the chairs. Expanding 
each term to order and recombining them we get 




+ a" 


— (cq + 20 ci + 4c2) {D^ + 2 DyD'j^ + (c 2 — ci)D^Di, 


+ C2 (SDpDjy — ADpDyDp + 2DiyDp^ 


G,p}+0{a^), (3.27) 


where the arguments (t, x) on the RHS have been omitted. Collecting ah results we dehne 
the expansion coefficients 


where the leading term dehnes the continuum limit, 

x) = dtBp{t, x) - (co + 8 ci + I 6 C 2 ) E BpGiip(t, x). 


(3.28) 


(3.29) 


Hence the correct normalization to reproduce the Yang-Mills gradient flow equation (1.1) 
is Co + 8 ci + 16 c 2 = 1, which we use to eliminate cq in the higher order terms: 



(3.30) 


+E 


(ci - C 2 )DID^ - C 2 X] (3-DpE)^ - iDpD^Dp + 2 D^D^p) 


G,p. (3.31) 


Before proceeding we remark on the presence of odd powers of a in the expansion, which 
seems at odds with our expectation that only even powers of a occur in this theory. The 
resolution of this apparent contradiction lies in the fact that the lattice helds Fp_{t,x) and 
Lp{t,x) should be dehned on the lattice link connecting x and x -)- a/i, rather than at the 
lattice site x. In Appendix B we demonstrate how the covariant re-expansion about the 
midpoint of the link, x = x + \ajl, eliminates such terms. While this problem will not 
affect our discussion of the 0 {a'^) counterterms, it clarihes that the corrections terms are 
indeed of order a‘^. 

We now proceed and work out the simplifations due to the held equations for x) 
and Lp{t, x). Varying the continuum action with respect to Lp one obtains the Yang-Mills 
how equation (1.1), whereas the variation with respect to Bp{t,x) yields 

dtLp = {^DpD,yLy + D^Lp) . (3.32) 
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Using the flow equation eliminates the 0(a) term and this is the reason why the 

O(a^) terms remain unaffected by the symmetrization about the midpoint x, once the field 
equations are taken into account. From the continuum ffow equation we derive 


di DyGpfj, — (SDpDy — ADpDjjDp + 2Dj^D^) (3.33) 

U z/,p 

This allows to rewrite the 0(a‘^) term as follows: 

^) = E| - - "2) + Dl) 

+ (^ + Cl - C2^ DpD^ - C2dtD^ (3.34) 


From the corresponding O(a^) flow action, 

5'2,fi[-B,L] =-2 j dt j d'^x^tr |L^(t,x)F^^)(t,x)| , 


(3.35) 


one may now directly read off the counterterm structures Qi that correspond with a given 
choice of the coefficients ci^ 2 - Unfortunately, there does not seem to be a choice such that 
52,11 vanishes. We also attempted to use Eq. (3.32) as follows: considering the term 


2C2 



ll,U 


(3.36) 


one may perform an integration by parts with respect to t. This generates a surface term 
at t = 0, 

- 2c2 / d‘^x E {Fp(t, x)Dj^Gi/^(t, x)}\^_Q , (3.37) 

which re-defines a coefficient of the counterterms entering 82 ^ (cf. Subsect. 3.5). Eq. (3.32) 
then leads to space-time derivatives acting on Lp, which can be integrated by parts (no 
surface terms are generated here) to redefine Fp . Unfortunately, this does not yield a 
solution with 52 ,h = 0 either. We notice, however, that 52 ,h with the Liischer-Weisz choice 
of coefficients ci = —1/12 and C 2 = 0, has a rather simple structure, 


*^2,01LW “ 



(3.38) 


To cancel this term is relatively straightforward. Starting from the lattice gradient force 
defined with the Liischer-Weisz action, 5 lW) we simply act with, 

1 + (3.39) 

on this gradient force, which yields the “Zeuthen flow” equation (1.5). The flow action 
52,11 for the Zeuthen flow does indeed vanish, i.e. we have successfully implemented 0 ( 0 ?) 
improvement in the 4-|-1-dimensional bulk. 
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3.4 0{a^) improvement of E{t,x) 


We here consider only the simplest observable, namely the action density E{t, x) of Eq. (1.4) 
The two most popular lattice discretisations of E{t, x) are referred to as plaquette (pi) and 
clover (cl) definitions, respectively. They are either obtained from the Wilson plaquette 
action or based on the so called clover leaf definition of the field strength tensor, 

~ x) + x) + R^y{t, x) + 3^))j4^JJ ) (3.40) 

which uses the 4 plaquettes (3.12), (3.17-3.19) in the fi — i' plane. The plaquette and clover 
lattice versions of E{t,x) are now given by: 


£'P'(t,x) = -^a \tT [P^v{t,x) + -2N 

E^\t,x) = -^'^tT{G‘^^^{t,x)G^^^{t,x)} . 


(3.41) 

(3.42) 


Pushing the classical a-expansion of the plaquette (3.12) to 0(a®) one obtains 
EP\t, x) = x) + ^ [^tr {D^G^i,{t, x)f + tr {DyG^y{t, x)f 

ll^U 

X] {G^y{t, X)f 

~ Y1 W + tr {Gf,u{t, x))^ + O(a^), (3.43) 


with the continuum limit P™“*(t,x) given by Eq.(1.4). Proceeding in this way for all 4 
plaquettes of the clover leaf we obtain the classical expansion 

E^\t, x) = E“°*(t, x) + ^ [^tr {D^G^y{t, x)f + tr {DyGf,u{t, x)f 

fl,U 

{dl + dl) tr (G^^(t, x)f + O(a^). (3.44) 

Several remarks are in order. First, the a-expansion of the plaquette yields contributions 
at every order in a, whereas the symmetries of the clover definition imply only even powers 
of a. The odd powers of a could be eliminated by averaging over the 4 plaquettes of the 
clover leaf, which, due to the trace operation, coincide with tr[P^;y(f, x)] for appropriately 
displaced arguments x. Second, note the total derivative terms which may appear at any 
order in a. Such terms do not contribute to the expectation value {E{t,x)), provided 
that the chosen set-up is translation invariant. This would e.g. be the case in a finite 
volume with periodic or twisted periodic boundary conditions, and thus in the limit of 
infinite volume. However, translation invariance no longer holds with either Dirichlet or 
Neumann conditions® as required for the Schrodinger functional [13] or with open boundary 

®Such boundary conditions are often imposed in the Euclidean time direction, combined with periodic 
boundary conditions in the spatial directions. In this case one may distinguish between the electric and 
magnetic components of E{t,x). In the latter, total derivatives only appear in the spatial directions and 
thus do not contribute to the expectation value. 
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conditions [14]. Similarly, when considering higher correlation functions such as the 2-point 
correlator of two fields E(t, x) total derivative terms cannot be ignored. We will here focus 
on the translation invariant case and from now on consider such total derivative terms 
negligible. This eliminates all the odd powers of a in the expansion of E'^^{t,x). Hence, 
both discretizations are on equal footing and counterterms O 2 for E^^ and E^^ are now easily 
identified as the O(a^) coefficients in the classical expansion. Given both a-expansions we 
observe that the O(a^) terms have the same structure, with the coefficients in the clover 
definition being larger by a factor of 4. In any case we observe that the linear combination 

x) = ^EP\t, x) - ^-E^\t, x ), (3.45) 

defines an O(a^) improved observable for which O 2 vanishes. An alternative O(a^) im¬ 
proved definition of E{t, x) can be obtained from the action density of a tree-level improved 
lattice action such as the Liischer-Weisz action (Eq.(2.1) with cq = 5/3, ci = —1/12 and 
C 2,3 = 0)- Here again, any ambiguity in the definition of a density from the action amounts 
to total derivative terms, which we consider negligible in the present context. 

3.5 Determination of 52,6 

In this subsection we list the gauge invariant local fields of dimension 6 which may appear 
in the boundary action 82 ^ of Symanzik’s effective action. Disregarding total derivative 
terms with respect to the space-time coordinates x, we find the following list of 7 candidate 
counterterms, 


Oi{x) 

= '^'^^{[D^Fpy{x)]Dp,Fpy{x)], 

(3.46) 

02ix) 

= ^ tr{[D^Fj.p(x)]D^F^p(x)}, 

(3.47) 

Osix) 

= '^^^{[D^,Fpy{x)]DpFpy{x)}, 

(3.48) 

Oa{x) 

= tr{Lp{0, x)DyFyp{x)}, 

(3.49) 

0^{x) 

= ^tr{L^(0,a:)L^(0,x)}, 

p 

(3.50) 

Oq{x) 

= dtti:{Gp.u(t, x)Gfj,u{t, x)}|t=o, 

(3.51) 

Oj{x) 

= ^tr{L^(t,x)5tH^(t,x)}|t=o, 

p 

(3.52) 


where denotes the field strength tensor of the fundamental gauge field. 

Again we apply the field equations. The Yang-Mills flow equation implies 

dtGfj,y{t,x) = [D^DpGpy — DyDpGpp] , (3.53) 

p 

so that, after taking into account the boundary condition Gpy\t=o = Fpy, we have 

Oq + 4O3 = total derivative, O7 = O4. (3.54) 
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This eliminates Oqj. The field equation (3.32) is not useful here. However, a third field 
equation can be derived by varying the action at t = 0 with respect to the fundamental 
gauge field Technically this is best done by discretising only the flow time in the 

4+1 dimensional continuum action and taking the limit of continuous flow time in the end. 
The resulting field equation is® 

^ ^ DyFy^{x) = -L^(0, x). (3.55) 

^ V 

This equation leads to the relations 

O 5 = -4^4, O 3 = -g^O^ (3.56) 

Hence one may also eliminate 03^5 in favour of O 4 . 

At this point it is useful to recall the situation in the standard 4-dimensional theory [19]. 
In fact there is a 1-parameter family of O(a^) improved actions, which, to tree-level, are 
parameterized by Xp as follows: 

co = 5/3 —24xp, ci = —1/12 + Xp, C 2 = Xp, C 3 = 0. (3.57) 

Expanding the action classically, the free parameter Xp is seen to multiply the counterterm 
O 3 . The counterterm O 3 is thus redundant for the improvement of standard observables. 
In principle one may thus tune the coefficients (3.57) to achieve O(a^) improvement of both 
standard and gradient flow observables. In practice however, these coefficients define the 
gauge action used in the Monte-Carlo simulation and the corresponding effective coefficient 
of O 3 should be regarded as fixed. One therefore needs to find an alternative way to achive 
improvement, and we chose to implement the counterterm O 4 (cf. Subsect. 3.6). 

Finally, we remark that the use of the field equation (3.55) in the counterterm basis 
holds for counterterm insertions only up to contact terms, namely whenever the countert¬ 
erm argument coincides with the location of some field in the correlation function under 
study. Such contact terms are thus absent for gradient flow observables localized at strictly 
positive flow times. However, we expect these relations to hold more generally, i.e. even 
if some fields in the correlation functions are defined at zero flow time. In this case we 
expect that the contact terms which make the difference are of the same form as the 0 {a‘^) 
counterterms to the fields in the correlation function and therefore just redefine these coun¬ 
terterm coefficients. This parallels the discussion in ref. [26] of on-shell 0{a) improvement 
in lattice QCD with Wilson quarks. 

3.6 Summary of Section 3 and some practical considerations 

Sect. 3 contains the main results of this paper and may appear rather technical. We 
therefore provide a short summary and comment on the practical implementation of the 
lattice counterterm O 4 . 

®While the continuum derivation may seem rather formal we note that a lattice version of this equation 
can be derived directly from the e-regularized 4-|-l-dimensional lattice action by a variation with respect to 
the link field U^{x), followed by the limit e —>■ 0. 
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There is a natural way of interpreting the gradient flow as a 4+1-diniensional local 
quantum field theory. The flow time t plays the role of the coordinate in the fifth dimension, 
which only takes on non-negative values {t > 0). The dynamics of the theory in the bulk 
(t > 0) is completely fixed by the deterministic flow equation. The classical nature of 
the theory for t > 0 allows to implement the Symanzik improvement programme in a 
rather simple way: all O(a^) cutoff effects produced by integrating the flow equation can 
be eliminated via a suitable discretization of the flow equation, which can be determined 
by the classical expansion to O(a^). Similar considerations allow to define discretized 
flow observables that are free of O(a^) lattice artefacts. The only remaining 0 ( 0 ^) effects 
are generated by the action at the boundary t = 0 , and are genuine quantum effects. 
They correspond to the usual O(a^) counterterms (3.46-3.48) in the 4-dimensional action 
affecting all lattice observables. 

To implement an O(a^) improved lattice action one first has to choose an O(a^) im¬ 
proved 4-dimensional lattice gauge action which amounts to choosing coefficients cq-s in 
Eq. (2.1) appropriately. It is well-known how O(a^) improvement can be implemented 
at tree-level, and also to order in the case of the pure gauge theory [20]. In addition 
one needs to incorporate a lattice version of O 4 such as to cancel the insertion of O 3 on 
observables without changing the coefficients co- 3 . 

To achieve this we remind the reader that the 4-|-1-dimensional set-up is used only 
for the theoretical analysis, whereas in practice one integrates the gradient flow equation 
numerically and evaluates any observable such as E(t, x) along the flow. It turns out that 
the insertion of O 4 can be realized by a change in the initial condition at t = 0 for the 
gradient flow equation. Since in this case and 5^(0, x) are not the same we need to 

fix the integration variables in the 4-|-1-dimensional field theory. We choose to integrate 
over the fundamental gauge field ^^(x) and the flow field i3^(t,x) for t > 0. Therefore on 
the lattice we choose to integrate over U^{x) and V)i(t,x) for t > 0. A shift in the initial 
condition can be implemented via 

V^it,x)\t=o = ^ (3 58) 

where Ch is the free improvement coefficient, and 5'g[t/] any 4-dimensional lattice action. In 
the 4-|-1-dimensional formulation with e-discretized flow time, the fields 1^(0, x) and U^{x) 
only enter in the terms 

Sg[Uf,] - x) 

where X^(t,x) is, up to terms of 0(e), the RHS of the flow equation. Now we can trade 
all references to V^{0,x) into Uf^{x), that is our path integral variable. Using Eq. (3.58) we 
can write 

U^(e,x)U;(0,x) = U^(e,x)C/t(x) - Cbgid.,^S^[U] + ... (3.60) 

where the dots represent higher order terms in the lattice spacing. Therefore the shift in 
the initial condition is equivalent (up to higher order corrections in a) to the insertion of 


-1 


fV),(e,x)V)J(0,x) - ij -eA^(0,x) | , (3.59) 
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the counterterm 


2cba^ ^ 04 ( 2 ;) = -2cba^ ^ tr {L^(0, x) {go^x,^,Sg[U]) } . (3.61) 

X X,fl 

Renaming the variable V)j(0,x) to gets us back to the previous situation with stan¬ 

dard boundary conditions, V)i(0,x) = U^{x) except for the extra O 4 term in the action. 
Hence we have successfully traded the modified boundary conditions for the flow equation 
for the O 4 term in the lattice action. In the next section we will determine its coefficient 
Cb at tree-level of perturbation theory. 

4 Perturbative analysis 

In this section we will study the Symanzik O(a^) improvement of the gradient flow in 
perturbation theory. This will allow us first to determine the improvement coefficient 
Cb to tree-level. Second, the study of the Zeuthen flow both in small volumes and for 
different observables will allow us to check explicitly that the use of a tree-level improved 
action for the simulation together with the tree level value' of Cb, the Zeuthen flow and 
a classically improved definition of the observable yields expectation values that are free 
of 0{a^) effects at tree level. As observables we choose first E{t,x). The contributions 
by the action, flow and observable to the cutoff effects of {E{t,x)) at tree level have been 
computed recently [22, 27]. Here we will show that the O(a^) tree level cutoff effects 
are absent not only in infinite volume, but also in a finite volume with twisted periodic 
boundary conditions, where the additional scale L leads to more stringent tests. Second 
we will consider the connected correlation function for E{t, x)E{s, y) and show that O(a^) 
improvement by the flow is also obtained in this case. 

4.1 Gauge fixing 

In perturbation theory one parametrizes the links in a neighbourhood of a classical config¬ 
uration as follows 


U;_i(x) = exp(agoAf,(x)); x) = exp{agoBij,{t, x)). (4.1) 

Note that this standard convention implies a re-scaling of the fields, 

y goAfj^, >■ qqB^^ ( 4 - 2 ) 

compared to the preceding sections. In perturbation theory it is convenient to use gauge 
symmetry to simplify explicit computations. In the context of the gradient flow, gauge 
fixing is performed by studying the generalized flow equation 

dtB^^\t,x) = D\^^G^^){t,x) + aD^^)d,B^J^\t,x), b\^\q,x) = A^{x). (4.3) 

^Note that all the improvement coefficients Ci and Cb have a perturbative expansion of the form c((?o) = 
-I- -I- . ... Since we are only concerned with tree level improvement we will omit the superscript 

in all improvement coefficients. 
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The superscript (a) serves as a reminder that covariant derivatives and field strength are 
made of the modified flow field Bj^\t,x), i.e. the solution of the above equation. Note 
that the original flow equation is recovered by setting a = 0. The key observation is that 
gauge invariant observables are independent of a [1, 3, 24]. In order to see this, one only 
has to check that the gauge transformation 

+ Ad^A -^, (4.4) 

where 

dtA = aAd^B ^,; A|^^q = 1, (4.5) 

transforms a solution of the flow equation with arbitrary a into one with a = 0. 

On the lattice the procedure is completely analogous. We consider the generalized flow 
equation 

a'^dtV^{t,x) = ^a'^\/J^\^\(^t,x)^dtA{t,x)]'^V^{t,x), (4.6) 

or, for the case of the Zeuthen flow, 

a^dtV^{t,x) = gl |- + a^V;^[A'l'(t,x)9tA(t,a;)]| V)f(t,x), 

(4.7) 

with initial condition with 1/^(0,x) = U^{x). One then easily verifies that the gauge 
transformation 

V^{t, x) = A(t, x)V^{t, x)A{t, X + /l)^ (4.8) 

transforms a solution with an arbitrary function A(t, x) into one with A = 1. A natural 
choice for the function A{t, x) then is given as the solution of the equation, 

A'^^tA = ad*^Bf_,{t, x ), A|^^q = 1. (4.9) 

Note that this is a particular application of the 4+1-dimensional gauge transformations 
described in Subsect. 2.3 and it is thus clear that gauge invariant observables remain unaf¬ 
fected by the choice of a. This can be turned around to provide checks on the correctness 
of a given calculation. In the following we drop the indices (a) (or A) from the fields and 
we will quote any intermediate results in Feynman gauge (a = 1). Some elements used for 
our checks of gauge parameter independence are given in Appendix C. 

4.2 Determination of Cb to tree level 

We first assume that the lattice is infinitely extended and expand the general class of 
actions, Eq. (2.1), to leading order in the coupling*^, 

Sglc'; {+’}] = 5 E / +{-p)+t’(p; A)+(J>) + O(9o), (4.10) 

IJ.,u ''P 

®See Appendix A for a summary of our notation and conventions. 
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where A is a gauge fixing parameter and explicit expressions for the lattice kernels, (p; A), 
are given in Appendix A. 

Similarly, the flow equation contains the gradient of a lattice action which, to leading 
order in the coupling, is parameterized by another action kernel, (p;a). The flow 
equation to this order then takes the form of the heat equation, 




(4.11) 


The initial condition for the flow equation Eq. (3.58) reads to leading order in the fields®. 


B^l{^,p) = ^ [V + a^CbKj^}{p-, 0) A^{p ), 


(4.12) 


(i) 

where is yet another action kernel. No gauge fixing term is required here, so that 
the gauge parameter is set to zero. The linearized flow equation (4.11) can now be solved 
easily 

B^^{t,p) = E H^.u{t,p■, a) [6,p + a^CbKl}}{p- 0)] i“(p), (4.13) 

v,p 

where is the heat kernel given by 


Hfj,i,{t,p-,a) = exp (^-tKB)(p-^a)^ 


{lU 


(4.14) 


Note that we have used here (p; a) as a matrix with respect to the Lorentz indices and 
the exponential has to be taken of that matrix. In the following we will often make use of 
such a matrix notation, in order to avoid an abundance of Lorentz indices. 

Finally, the observable E{t,x), being an action density, can be parameterized by a 
further lattice action kernel, K^°\p,0), with gauge fixing parameter set to zero. To this 
order we then obtain for the expectation value 

{E{t,x)) = ^ ^ J Tr|A:(°)(p;0)Z)(t,t,p;a,A) I + C>(po), (4.15) 

where the trace is over Lorentz indices only and the gauge field propagator at positive flow 
time is defined by 


{Ki^,p)Bl{t,q)) = {2Tx f 5^^\p + q)5‘^'’D^^{p,s,t\ a, X). 


(4.16) 


Due to the relation (4.13), this propagator depends implicitly on both gauge parameters, 
a and A, of the flow equation and of the the action, respectively. Introducing the standard 
4-dimensional gauge field propagator 


{A;{p)Ai{q)) = (27r)45(4)(p + q)5-^D^,{p- A), 
this propagator is the matrix inverse of the action kernel, 

a(“)(p,A)D(p,A) = 1, 


(4.17) 


(4.18) 


®Note that higher orders in the helds imply higher order in the coupling, too, cf. Eq. (4.1). 
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and the gauge fixing parameter A must be non-zero for the inverse to exist. Using these 
ingredients, the gauge field propagator at positive flow time can now be written as follows, 


D{p, s,t] a, X) =H{s,p;a) t + a^ci,K^'^\p;0) 


1 T 


x£)(“)(PjA) 11-|-p; 0) H{t,—p-,a) 


(4.19) 


where we have denoted the matrix transpose by the superscript T. 

In summary, the choices of action, flow and observable discretization correspond to 
the choice of three action kernels. Finally the shift in the initial condition is encoded in a 
fourth choice of kernel. Explicit expressions for some popular choices of kernels are given 
in the Appendix C. 

In order to obtain the leading order cutoff effects we now expand the kernels as follows. 


K{p; A) = A) + a^Rip] A) + O(o^), (4.20) 


where the continuum kernel is given by 

V - (1 - ■ (4-21) 


Using the continuum kernel only and neglecting cutoff effects we thus obtain the well-known 
continuum result in infinite volume, 

(E(t, X)) = + 0(9*. a\ £„““(*) = ^[28 Jv - 

Explicit expressions for the correction terms R^y{p;X) are given in Appendix C. In order 
to compute the leading correction to the propagator D^^{p-,X) and to the heat kernel 
H/iiyit,P',a) it is convenient s to work in Feynman gauge (A = a = 1), since in this case 
K'^^{p] 1) is proportional to df^u- Working in a general gauge is however not much more 
difficult and serves as a check that the gauge dependence actually cancels in the final 
evaluation of the observable. A few technical details pertaining to such a check are given 
in Appendix C. 

In the following we will use Feynman gauge and remove the gauge parameters as 
arguments of the action and flow kernels. We will also omit them in the kernels for the 
observable and initial conditions however, with the understanding that they must be set 
to zero in these cases. In Feynman gauge (A = a = 1) it is straightforward to check that 


Dfxuip) 


H^xu{t,p) 


p^ 


+ 0(a^ 


e - a^tRfj,uip)] + 0(0“^) 


(4.23a) 

(4.23b) 


and finally, putting all the pieces together and after some algebra, we get 


£o{t)=£r\t){l + 


t - 


( 4 °^ - d^ 4 )'h ,-2 + (4°^ - 4“^ + 2Cb) J2,0 - 
- 2 d 4 \h,o - 24^^-^2,2] } + 0 {a^) , 


(4.24) 
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Discretization 

di 

d2 

Plaquette 

-1/12 

0 

Liischer-Weisz 

-1/36 

1/36 

1 Plaquette — | Clover 

-1/36 

1/36 

One-parameter tree-level improved 

-1/36 

1 /36 — Xp 

Clover 

-1/4 

-1/12 

Zeuthen 

0 

0 


Table 1: Values of the coefficients in the O(a^) terms of t‘^{E{t, x)) in infinite volume. The 
one-parameter family of tree-level improved actions corresponds to the choice of coefficients 
Eq. (3.57), the Liischer-Weisz tree level improved action being the particular choice with 

Xp = 0 . 


where the constants Jn m are defined by 


T —An+m)/2->P^ ' 

^ f g-2tp2 


and 


EjPMr ^>0 


p = 


1 -1 


EniPp)"" ri<0 


In fact it is straightforward to evaluate the integrals with the result, 


(4.25) 


(4.26) 


<74,-2 — 1/2, <72,0 — 1) <74,0 — 3/4, J2,2 —3/2. (4.27) 

The coefficients must be independent of the gauge parameters a and A and we have 

checked this explicitly. Their values depend on the choices made for the various kernels. 
For example, for a general action of the form Eq. (2.1) we have 

12 2 2 

- gCi -k -C 2 -k -C3 , (4.28a) 

d2 = • (4.28b) 

Table 1 summarizes the values of the coefficients loi’ most common choices. It is 

easy to see that the use of the Zeuthen flow together with the tree-level improved Liischer- 
Weisz action and any classically improved discretization for the observable (see section 3.4) 
has no tree-level 0{a?‘) cutoff effects as long as = 0. Therefore, to tree-level, the Liischer- 
Weisz action (ci = — 1/12,C2 = 0) produces tree-level improved results for gradient flow 
observables. For the case of a generalized tree-level improved action Eq. (3.57) we have to 
choose ^ 

Cb = -2^P’ (^-29) 

in order to obtain tree-level improvement. 
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As the reader can see, besides the Zeuthen flow there seem to be many ways to cancel 
the tree-level O(a^) effects (see also [22]), as these are encoded in a single term, once the 
numerical values (4.27) and for d \^2 (cf- Table 1) are inserted into Eq. (4.24). We are 
thus led to look for more stringent tests of 0 {a^) improvement by looking at a variety of 
observables and/or kinematics. After all, rather than improving a particular observable in 
a specific situation (e.g. in infinite volume), Symanzik improvement is designed to work for 
any observable in both finite and infinite volume. 

4.3 Twisted periodic boundary conditions 

A stringent test of our computations can be made when studying t‘^{E{t,x)) in a finite 
volume. Due to the presence of a new scale L, the cutoff effects will in general depend on 
the dimensionless ratio c = y/^/L. Improvement requires that the tree-level cutoff effects 
vanish for all values of c = y/M/L. 

As a finite volume renormalization scheme, we will use twisted boundary conditions 
for our gauge field. In this setup, the gauge field changes by a gauge transformation when 
displaced by a period. Gauge invariant quantities are still periodic, but the absence of 
zero-modes in the perturbative expansion turn out to be very convenient for our analytic 
computations. The gradient flow has already been studied in this setup, and we will not 
give much details here but refer the interested reader to the work [15] and the references 
cited therein. 

We will only need the perturbative expression of {E{t,x)) to leading order, given by 


{E{t, x)) = go£o{t, c) + 0{gQ) (4.30a) 

with 


[l + a^c,K(^\P)] D^-\P) 

p 

1 + (-P) ^ (t, -PfK^°'^ (P) I . (4.30b) 


Note that the expression is almost identical to the infinite volume one, except that the 
momentum integral has been substituted by a sum (hardly a surprise). The particularities 
of the twisted boundary conditions are hidden in the sum and momentum symbols. First 
notice that the momentum (with capital letters P^) can be uniquely decomposed as 




27rnp 

L 


+ 


27rn^ 
NL ’ 


with = 0,..., L/a — 1 and 


= 


0, if /U = 0, 3 , 

0 ,...,iV-l, ifM = l,2, 


(4.31) 


(4.32) 


i.e. there is the usual space-momentum, but in the directions of the twisted plane xi — X 2 
the momentum lives in an apparently larger lattice of size NL, where N is the rank 
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of the gauge group. Finally the sum symbol ^ means sum both over and h^, but 
without the terms with hi = h 2 = 0. In particular the sum has no term with a zero total 
momentum. Notice that the colour factor N'^ — 1 is produced by the sum over h^. 

The algebra is very similar to the one of the previous section, with the important 
difference that now the sums actually depend on the dimensionless ratio c = y^/L. In 
fact fixing the flow time in units of the volume in this way we get 


,2 , 


£o{t, c) = £r\t, C){1 + ^ {d[°^ - 4“V4,-2(c) + - dr + 2cfe) J 2 ,o(c) - 


(o) .(a) 


L 


-2^^^ J4,o(c) - 2drJ2,2{c) } + 0{a^), 


Af) 


(4.33) 


where 

= Y^^i(0|^7rc2) [hi (OlzTTcViV^) - 

and the third Jacobi theta function reads 


Finally the functions J7ij(c) are 


function reads 
hg {z\t) = 

n. 

given by^^ 


> ( 4 - 34 ) 


(4.35) 


•Ji,j (c) 


f CTT 




\*+i ^ exp{—c 27 r 2 (n + h/A^)2] 


(4.36) 


In the limit c = y/^/L —)• 0, we recover the expressions of the infinite volume, in 
particular 

lim J7)j(c) = Jij , (4.37) 

c^O 

but for non-zero c the functions J7ij(c) are in general linearly independent (see Fig. 2). The 
coefficients still the same, and the reader can check that the tree-level 0 (a‘^) 

cutoff effects given by expression Eq. (4.33) vanish for all values of c when one uses our 
improved setup (i.e. Liischer-Weisz action, Zeuthen flow and Liischer-Weisz observable). 
Any other choice of improved action toghether with the appropriate choice of Cf, also does 
the work. For this to happen it is crucial that the flow coefficients d) 2 both zero, since 
the functions J4,fi{c) and J 2 , 2 {c) are linearly independent. In particular it is easy now to 
check that the so called Symanzik flow in the literature [28] or any set of coefficients in [22] , 
does not remove the tree-level cutoff effects in finite volume. For the Zeuthen flow both 
coefficients do identically vanish, so that O(a^) effects are indeed removed as expected on 
theoretical grounds. 


4.4 The connected 2-point function of E{t,x) 

Further tests of the Zeuthen flow can be obtained by considering different operators at 
positive flow time. In particular, we now consider the 2-point function of E{t,x) in a 

^^Negative powers (j < 0) have to be understood as in the infinite volume, Eq. (4.26). 
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1.2 



Figure 2: Ratio of the finite volume functions Ji^j{c) (Eq. (4.27)) over the infinite volume 
predictions Jjj (Eq. (4.36)). When c > 0.2 there are significant differences between them. 
Moreover the different functions Jij{c) are in general linearly independent. 


periodic box of size L with twisted periodic boundary conditions (i.e. the same setup as 
above), 

t^s^{E{t, x)E{s, 0))c = [{E{t, x)E{s, 0)) - {E{t, x)){E{s, 0))] . (4.38) 

The factor renders this quantity dimensionless, so that it can be considered a function 
of the two dimensionless parameters, 

(4.39) 

1j 

Computing to leading order in the coupling the result can be written in the form, 

t^s^{E{t, x)E{s, 0))c = t^s^glM{t, s; x) + 0{gl) , (4.40) 



with, 

s-,x) = ^^ Y. Trj K{P, Q)H{P, t)D{P)H{-P, sf 

xK{Q,P)H{-Q,tfD{QfH{Q,s)^ . (4.41) 

The generalized kernel K{P,Q) encodes the discretization of the observable. Up to terms 
of O(a^) it is given by 

K^,{P, Q) = Q) + O(a^), (4.42) 
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with the continuum kernel given by 


Q) = E V - • (4-43) 

p 


The finite volume calculation for {E(t,x)) has taught us that the 0(a‘^) contributions of 
the flow have to cancel by themselves, i.e. a cancellation with other O(a^) contributions 
by the action or the observable are not possible, due to the linear independence of the 
momentum sums. In order to assess the improvement of the Zeuthen flow it is therefore 
enough to focus on these O(o^) contributions. Using again the Feynman gauge for flow 
and action, we obtain 0 ( 0 ^) terms from the flow of the form. 


2 ^{P+Q)x -(t+s){P^+Q^) 

1024 ^ 


P,Q 




X Tv\ Q) {tR{P) + sR{Pf) P) \ , 


(4.44) 


and a second term with a similar structure. In both cases it is useful to note the property 
of the kernel, 

Q) = T{Q)K^°'^\P, Q)T{P), (4.45) 

where T{P) is the transverse projector, 

V(P) = V - ^ ■ (4-46) 

The 0{a‘^) correction to the Zeuthen flow kernel, R^{P), has the nice property that 

T{P)R^{P)T{P) = 0. (4.47) 

Hence we can conclude that the Zeuthen flow does not contribute any 0{a'^) effects to this 
2-point function either. Due to the different Lorentz index structure of this case compared 
to the simpler case of E{t,x), and to the fact that now, in general, the cutoff effects 
are functions of two variables {c,d), this test imposes further constraints on the possible 
improvement solutions. In particular, the so called chair flow in [27], which happens to 
also cancel the 0{a‘^) effects of {E{t, x)) in a finite volume, can be shown to produce 0{a‘^) 
contributions to the 2-point function considered here. 


5 Conclusions and Outlook 

We have systematically investigated the structure of 0{a‘^) effects in flow quantities using 
Symanzik’s approach applied to the 4-|-1-dimensional local formulation of the theory. Im¬ 
provement to 0{a‘^) for gradient flow quantities appears to be easier than one might have 
thought, mainly due to the classical nature of the gradient flow equation. In particular 
the classical a-expansion is sufficient to obtain the counterterms for both local composite 
operators at positive flow time and the action in the 4-|-1-dimensional bulk (i.e. due to the 
absence of loops in the bulk, no new counterterms are generated). 
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Our main results are summarized in the Zeuthen flow equation (1.5) and the improved 
lattice definitions of the observable E{t,x), either as linear combination of clover and 
plaquette dehnitions (3.45) or as the action density of the tree-level improved Liischer- 
Weisz action. We have shown that the integration of this Zeuthen flow equation and the 
evaluation of classically improved observables do not produce any O(a^) effects to any order 
in the coupling or, indeed, non-perturbatively. At this point it is important to remark that 
although the analysis has been performed in the context of pure gauge theories, due to 
the classical nature of the flow equation, the aforementioned results are still valid in QCD 
or if any number of fermions in any representations are coupled to our gauge field. In 
the particular case of the pure gauge theory the only C>(o^) effects originate either from 
the 4-dimensional lattice action or from the additional counterterm parameterized by Cb in 
the modified initial condition (3.58). Tree-level O(a^) improvement is achieved with the 
Liischer-Weisz gauge action and Cb = 0- 

We have explicitly checked that the proposed Zeuthen flow equation does not gen¬ 
erate any O(a^) contribution to tree level for a variety of gluonic observables (different 
observables in arbitrary volumes). In doing so, we have shown that other proposals of the 
literature to improve the gradient flow (i.e. the r-shift in [21], the coefficients in [22] or the 
chair flow in [27]) in fact do produce 0{a‘^) effects in some of the considered observables. In 
this sense, these proposals only produce vanishing O(o^) cutoff effects in some particular 
situations (i.e. {E{t,x)) in inhnite volume), and this cancellation should be regarded as 
accidental, and not as improvement. 

Our results can be extended in various directions. First, it appears straightforward 
to extend the classical a-expansion to further observables, for example the the energy- 
momentum tensor. When considering n-point correlation functions of such observables with 
n > 1 or if boundary conditions do not respect translation invariance in some directions 
(as is the case with SF and open boundary conditions), some additional work is required 
to also eliminate total derivative terms which may contribute at any order in a. We also 
note that the improvement of observables and the flow equation are conceptually separate 
from the O(a^) effects at t = 0. It is therefore conceivable to push the expansion further, 
in order to also cancel terms at O(a^). It is not clear how complicated this would be for 
the flow equation, but it is certainly an option for observables. However, one should be 
aware that higher order improvement would typically render these observables less local in 
lattice units. Another natural generalization would be the extension of our work to include 
fermions and the fermionic flow equation, introduced in ref. [4]. 
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A Conventions and notation 


We will use the summation convention for colour indices 


a,6,... = 1,..., - 1, 


(A.l) 


but not for space-time indices as this may lead to confusion in the discussion of 

lattice artefacts. Trace over color indices will be denoted by tr (lower case), while trace 
over Lorentz indices will be denoted with the symbol Tr (upper case). 

SU(N) gauge fields live in the Lie algebra su(A) and are traceless antihermitian A x A 
matrices. Any element X £ su(A) of this algebra can be written as A = A“r“ where the 
components A“ are real numbers and the generators r“ are themselves antihermitian A x A 
matrices chosen to obey the normalization 

tT{T^T^) = -^6ab. (A.2) 

On the lattice the links U^{x) belong to the gauge group SU{N). For an arbitrary 
function of the link variables f{U^{x)), the Lie-algebra valued derivative is given by 

= T-dlJ{U^{x)) = r“A/(e^T“[/^(^)) . 

e=0 

Fourier transformations on an infinite lattice with lattice spacing a are defined as 

A^{x) = [ 

Jp 

where 

f _ d^p 

Jp J-n/a (2vr)4 ■ 

On a hypercubic lattice of volume we define 

= , (A.5) 


(A.3) 

(A.4) 


with Pf^ = 27rn^/L and = 0, L/a — 1. It is convenient to introduce the lattice 
derivatives 


(/)(x afi) - (pix) 

d^4>{x) = -, 

a 

^ </>(x) - (/)(x - a/i) 

= - , 

^ a 


(A.6) 

(A.7) 
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and also the covariant derivatives given by 


= U^{x)f{x + ail)U^{x)^ - f{x) , (A.8) 

= f{x) - U^{x - afi)^f{x - afi)U^{x - afi). (A.9) 

(A.IO) 

Along the work we use the following definitions of lattice momenta 

2 

= -s\n{apf,/2), (A.ll) 

= - sin(ap^), (A. 12) 

a 
1 

= -cos(ap^/2). (A.13) 

a 


B Absence of odd powers of a in the classical expansion 

In this appendix we demonstrate that the apparent presence of odd powers of a in the 
classical expansion (3.28) is an artefact of the way the expansion was set up. In particular, 
we will show that Symanzik’s effective action for the flow action only contains terms which 
are even powers of a. 

B.l Re-exanding around the midpoint of the link 

Indeed, the expansion about x does not account for the fact that the equation is derived 
for a given link variable Vfj,(t,x), relating the lattice points x and x + afi. Odd powers 
of a in the expansion are due to this asymmetric treatment, as can be shown explicity to 
all orders in a for the LHS of the flow equation, Eq. (3.22). First, we define the unitary 
matrices n^(t, x) as the parallel transporters along the half link from x to the midpoint 

X = X + ^ajl, 

0^(t, x) = "P exp |a ^ dtt (t, z(tt)) I (B.l) 

i.e. compared to the path ordered exponential V^{t,x) (3.8) we here only integrate over 
second half of the path parameterizing the link. Now we can perform the parallel transport 
to the midpoint x, defining 

Ff,{t,x) = n^(t,x)F^(t,x)Q^(t,x)“\ (B.2) 

(B.3) 

and, analogously, L{t,x), such that the term in the lattice flow action density, 

= tr{L^(t, x)F^{t, x)} = tr{Z^(t, x)F^{t, x)}, (B.4) 

can be expressed in terms of fields defined at the midpoint. To obtain the expansion in 
powers of a about x one may simply re-expand the expansion about x obtained previously. 
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The parallel transporters x) then merely render the Taylor expansion covariant. Pro¬ 
ceeding in this way yields, for the first term of 

= dtB^,{t,x), (B.5) 

n=0 1 ^ >' 

which explicitly contains even powers of a only. For the gradient force term in F^(t,x) 
we have only worked out the hrst few orders of the a-expansion explicitly. Therefore, the 
re-expansion cannot be carried out to all orders in a and it is thus advisable to resort to 
some more general argument based on symmetries. 


B.2 Reflection symmetries 

We now consider the flow action in Eq. (3.23), but restricted to plaquette and rectangle 
terms, as this is sufficient to discuss the case of the Zeuthen flow. We consider a coordinate 
reflection TZa in direction a. The point with coordinates transforms into x'^ with 


. x^ 

The gauge field transforms under TZq 

Vf,{t,x) — 


= 


—Xo, if /r = a. 


(B.6a) 


Va{t, x' — aa)^, if /i = a, 

y^{t,x'), if/iT^a, 


(B.6b) 


One may then show that the gradient force terms (the RHS of the flow equation), transform 
for plaquette and rectangle terms, as follows: 


x) —> 


—Va{t, x' — ay^ait, x' 

Xy,x') 


aa)Va{t, x' 


aa) if ^ = a, 
if ^ 7 ^ a. 


(B.6c) 


In fact, the same transformation behaviour is found for the left hand side of the flow 
equation, so that Eq. (B.6c) equally holds for Fy,x) of Eq. (2.5). Hence, if the same 
transformation behaviour (B.6c) is imposed on the Lagrange multiplier field, L^(f, x), we 
obtain for the different parts of the action density. 


tr (L^(t,x)F^(Lx)) 


tr {La{t, X — aa)Fa{t, x' 
tr {Ly,x')Fy,x')) 


aa)) if // = a, 
if /X 7 ^ a. 


(B.7) 


In particular, the action is invariant under such a reflection, as the only effect consists in 
a re-ordering of the terms in the sum over the Xo-coordinate^^. 

^^One may think of the inhnite lattice as a limiting case of hnite lattices with periodic boundary conditions 
where the sum over Xa is finite and the re-ordering of terms in the sum is unproblematic. 
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B.3 Example: reflection TZa of tho Wilson gradient force 

It is instructive to derive Eq. (B.6c) for the case of the plaquette action in some detail. 
The gradient force in this case has the form, 

x) = ^ x) + (B.8) 

and we need to distinguish the two cases ^ = a and /r / a. Starting with /r = a and 
setting y = x' — aa WQ obtain the transformation behaviour of these plaquettes 

Pav{t, x) — Va{t, y)^Pau{t, y) V«(t, y), (B.9a) 

and, similarly, 

Qau{t,x)^ —> Vait,y)^Qau{t,y)Va{t,y). (B.9b) 

Summing both expressions and taking the antihermitian part we thus obtain 

x) -|- Qau Pau{ti y) + Qah'1 Pa{t, J/) , (B.IO) 

where we have used the relation (M1 )ah = — valid for any square matrix M. 

Next we consider the case y ^ a. The transformations of the plaquettes in this case 

read 



^ iQ^ia{t,xy, 
\p^u(t,x'), 

if I' = a, 

if u ^ a, 

(B.lla) 

Q{f'l x)^ 

^ (Pf,a{t,x'), 
\Qf,u{t,x')^, 

if u = a, 

if V ^ a. 

(B.llb) 


Hence, Eq. (B.6c) follows and the part of the lattice flow action containing the Wilson 
gradient force is indeed invariant under a reflection TZa- We have also verihed that this 
remains true for any gradient force obtained from lattice actions containing both plaquettes 
and rectangles, such as the Liischer-Weisz action. 

B.4 Lattice vs. continuum reflections 

We now consider a total reflection, TZ = TZqTZiTZ2TZ3 of all space-time coordinates, i.e. 

TZ : X — )■ X = —X (B.12) 

The part of the flow action density for fixed index y then transforms as follows: 

ti{L^{t,x)Ff^{t,x)} —^ tr x' - ay)} . (B.13) 

It is not difficult to see that the 0(a) offset in this transformation is again an artefact of 
the asymmetric treatment of the links. In fact, deflning again the midpoint x = x + ^afi, 
and using the transformation of the ’’half link variables”, H(t,x) (B.l), 

TZ: H^(t,x )—> Qfj,{t,x' - (B.14) 
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we find that the transformation behaviour of the fields at the midpoint is given by 

Lf,{t,x) —^ -Lf,{t,-x), (B.15) 

F^{t,x) —^ -F^{t,-x), (B.16) 

i.e. the reflection TZ, once expressed in terms of the fields at the midpoint x takes the same 
form as its continuum counterpart. Therefore, the a-expansion of the corresponding part 
of the lattice flow action cannot generate terms that are odd under TZ, i.e. any term 

with an odd number n of Lorentz indices can be excluded. This together with the obser¬ 
vation that all Lorentz vectors L^, dfL ^,...) have odd canonical dimension, implies 
that any term containing an even number of them must be even dimensional and thus be 
accompanied by an even power of a. 

Treating the parts of the lattice flow action density with other values of fx in the same 
way, no odd powers of a can be generated in the expansions about the respective midpoints 
of the links relating x + afi and x. As these midpoints coalesce to a single point in the 
continuum limit this establishes this property for the a-expansion of the complete lattice 
flow action, for gradient force terms containing plaquette and rectangle terms. 

For these considerations to extend to the Zeuthen flow we only need to check that the 
correction term, 

v;v^x^it,x), (B.18) 

transforms like X^{t,x) itself under the reflection TZ. This is indeed the case, so that 
the absence of odd powers of a is confirmed for the Zeuthen flow, too. Finally, while it 
is plausible that these considerations extend to gradient force terms derived from lattice 
gauge actions containing the “chairs” and “parallelograms”, we did not check this explicitly, 
as it is not needed for the discussion of the O(a^) improved Zeuthen flow. 

C Action and heat kernels to O(a^) 

C.l Free lattice actions and their kernels 

The choice of observable, action and flow at tree level can be parameterized by the kernels 
of free lattice actions, i.e. the gauge action expanded to second order in the gluon fields, 
possibly supplemented by a gauge fixing term. If a generic lattice action with Wilson 
loops of length 4 and 6 is chosen then these are parameterized by a set of coefficients Cj, 
f = 0,1, 2, 3. An alternative is provided by directly inserting the clover leaf deflation of the 
gluon held strength tensor gives into a continuum like action density. In momentum space 
any of these actions is written 

Sg[U-,{c,}] = ^A“(-p)A^,(p;A;{Q})A“(p) + 0(ffo) (C.l) 

Note that for the case of a finite volume with twisted boundary conditions, the expressions 
of the kernels A) are unchanged, but the integrals over momenta have to be substi¬ 

tuted by sums and the momentum has to be interpreted as the sum of the space and color 
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Discretization 

Kpiuip'i X) 

Plaquette 

Liischer-Weisz 

Clover 

Zeuthen 

p‘^Sfj,y - {1- X)p^py 

P^5^y - (1 - X)p^,Py + [{p^ +P^pI)6^V - PtMPu(.Pl+pl)\ 

p'^C^S^U - PflCf^PyCy 

(1 - a^pl/l2) \^p^S^y - Pf,py -b [{p'^ + p^pl)6f,y - Pi^pyiPl + pD] 

1 -b Xp^py 

Discretization 


Plaquette 

Liischer-Weisz 

Clover 

Zeuthen 

- (1 - >^)lp^iPu{pl + pD] 
hp^pl^iiu - ^p,,pApI+pI) 

- ( 5 / + \p‘^pD + MP^iPu{pl+pl) 

jiP^Pu [(1 + x)pI - (1 - x)pI] 


Table 2: Kernels K^y corresponding to different choices of discretization, and discretiza¬ 
tion effect corrections for some of the most popular choices. See appendix A for any 
unexplained notation. 


momentums (see Eq. (4.32) and the subsequent discussion). Gauge fixing is performed in 
any kernel by adding the usual gauge fixing term 


K^y{p] A) = Kfj,y{p; 0) Xpf^pu . 

Expanding the kernels to O(a^) around their common continuum limit, 
K{p- A) = A) + a^R{p, A) + O(o^) 


(C.2) 

(C.3) 


with 

V + {^- '^)Pt^Pp (C.4) 

the leading cutoff effects are encoded in the structure of R^y{p;X). For example for a 
generic action made of an arbitrary linear combination of loops of 4 and 6 links Eq. (2.1) 
we have 


K^^\p) = p^5f,y - (1 - X)pf,py - a^{ci -C2- cs) [{p^ + p^pl)S^y 
- a^{c2 + Cs) [{p^fS^y - p^Pf,Py] , 


PtiPuiPl+Pl)] 

(C.5) 


and 


R^^'> = - 


+ P^^Pu 


^ + Cl - C2 - Cs ) / + (ci - C2 - C 3 )p'^pI + (C2 -|- Cs)(p^)^ 


-^flU 


1-A 

24 


Cl - C2 - Cs ) (p^ + pI) -b (C2 + C3)p^ 


(C.6) 


The expressions for K^y and Rfj,y for the other common choices of discretizations are 
written in table 2. Note that the Zeuthen flow equation, even if it is not derived from 
the gradient of an action, can also be parametrized to tree-level by a kernel. The main 
difference is that the property 

K^yip) = Ky^ip) , (C.7) 


that is obeyed by any kernel derived from the gradient of an action (i.e. a consequence of 
the action being real) does not hold for the Zeuthen flow. 
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C.2 Heat kernels and propagators to O(a^) 

Given a kernel K with arbitrary value of the gauge parameter A (or a in case of the flow 
kernel), we would like to obtain the correction term to either the propagator, i.e. the 
inverse of K, 

Dip-,\) = K{p-,X)-\ (C.8) 

or the heat kernel 

H{t,p]a)=exp{—tK{p]a)), (C.9) 

given the expansion of the kernel Eq. (C.3). 

Starting with the propagator, we can formally invert, 

D = + a^R + 0{a^)] “^ = [l + + 0{a^)] (C.IO) 

and then expand in to obtain 

D{p- A) = A) - a^D^°-^\p- X)R{p- A)D““*(p; A) + 0{a^). (C.ll) 

To work out the 0{a?) piece of the heat kernel we define the transverse and longitudinal 
projectors 

L,,{p) = (C.12) 

t- pZ f- pZ 

in terms of which 

a) = [T^,y{p) + aL^^{p)] , (C.13) 

and the heat kernel in the continuum is given by 

H^°’^\t,p-, a) = exp a)) = e-^P"T{p) + (C.14) 

Note that K^°'^^(p;a) and R{p]a) do in general not commute. Nevertheless, it is not 
difficult to work out the expansion to 0{a?‘) (i.e. to first order in R{p;a)). Inserting (C.3) 
in the exponent one obtains 

+ 0(^2), u = {l-a)ph, v =-ah. (C.15) 

Then, noting, that for n > 1, 

{uL + vR)^ = u^L + u^-^v {LR + (n - 2)LRL + RL} + 0{v^) , (C.16) 


the exponential series can be resummed with the result 


uL+vR _ sr^ (uT + vR)"' _ 


E 

n=0 


= T + e“L + u 


n! 


e“ - 1 / ^ 


R — R-\ -(i? + uLRL'j 


u 


+ 0(u2), (C.17) 


where 


R = LR + RL - 2LRL. 


(C.18) 
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The result for the heat kernel then is: 


H{t,p;a) = + a^te —i? 

1 — _ 'I 

+ - (1 _ [^ + (1 - a)tp^LRL] I + O(o^), (C.19) 

where we have left out the arguments for the sake of readibility. Note that the choice of 
Feynman gauge a = 1 for the heat kernel is not a problem, as the apparent singularity 
cancels, with the result 

H{t,p, ; 1) = {1 - a^tR{p-, 1)} + O(a^). (C.20) 
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